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Abstract 

For lnAs/GaAs(001) quantum dot (QD) system, the wetting layer (WL) evolution and its temperature dependence 
were studied using reflectance difference spectroscopy and were analyzed with a rate equation model. WL 
thicknesses showed a monotonic increase at relatively low growth temperatures but showed an initial increase and 
then decrease at higher temperatures, which were unexpected from a thermodynamic understanding. By adopting 
a rate equation model, the temperature dependence of QD formation rate was assigned as the origin of different 
WL evolutions. A brief discussion on the indium desorption was given. Those results gave hints of the kinetic 
aspects of QD self-assembly. 
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Background 

Epitaxial semiconductor quantum dots (QDs) have 
attracted much attention because of their application po- 
tential in novel optoelectronic devices [1-3]. They are 
usually fabricated utilizing the lattice mismatch between 
the epitaxial layer and substrate or the Stranski- 
Krastanov (SK) growth mode. It can be described as 
follows: for small coverage, two-dimensional (2D) layer- 
by-layer growth and pseudomorphic formation of wet- 
ting layer (WL) take place. When the WL reaches a cer- 
tain critical thickness (CT), a 2D to three-dimensional 
(3D) transition starts, and QDs form on the substrate. 
QDs with high homogeneity in their size and shape are 
highly advantageous in applications. Basically, the WL 
configuration would also influence the optical properties 
of QDs and the performance of QD-based devices [4-7]. 
A controllable growth of QDs with desired properties 
requires a comprehensive understanding on the growth 
process. Therefore, it is necessary to have a clear under- 
standing of the WL evolution during the QD self- 
assembly. 
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The commonly accepted thermodynamic understand- 
ing of the SK mode describes the QD formation on top 
of a WL of a certain thickness. However, it is not accur- 
ate in real situations. It has been reported that in the 
Ge/Si QD system, the WL thickness decreases after QD 
formation [8-10]. It is interpreted in the regime of kinet- 
ically controlled QD formation and growth. Since mater- 
ial transfer from WL to QDs sustains the QD formation 
and growth, a large material consumption rate by QD 
formation may induce the observed WL erosion [8,9]. 
As for InAs/GaAs system, a step erosion of WL has also 
been observed after QD formation [11,12]. Until now, 
there is no complete description of the WL evolution 
and its growth condition dependence. In our previous 
work, reflectance difference spectroscopy (RDS) was 
used to study the WLs in self-assembled nanostructures. 
Due to its sensitivity, heavy hole (HH)- and light hole 
(LH) -related transition energies before and after QD for- 
mation can be directly obtained from the resonant struc- 
tures in the spectra [13-16]. In this paper, we studied the 
WL evolution and its temperature dependence based on 
RDS measurements. We found that, generally, there 
were two kinds of WL evolution, with deposition de- 
pending on growth temperatures. They were well 
explained in the regime of temperature dependence of 



© 2012 Zhang et al.; licensee Springer. This is an Open Access article distributed under the terms of the Creative Commons 
Attribution License (http://creativecommons.Org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction 
in any medium, provided the original work is properly cited. 



Zhang et al. Nanoscale Research Letters 2012, 7:600 
http://www.nanoscalereslett.eom/content/7/1/600 



Page 2 of 6 



QD growth rate with a rate equation model The 
concave-up style of evolution was considered as a clear 
evidence for a non-zero QD growth rate when the WL 
thickness was smaller than the CT. We also gave a sim- 
ple discussion on indium desorption during self-assem- 
bly. All of these results showed the kinetic aspects of 
WL evolution in the SK growth. 

Methods 

Six InAs/GaAs(001) QD samples with different growth 
temperatures (from 490°C to 540°C, with an increment 
of 10°C) were grown in our Riber-32p molecular beam 
epitaxy (MBE) system. A 100-nm GaAs buffer layer was 
firstly deposited on 2-in semi-insulating GaAs substrates 
at 600°C. A nominal InAs amount of 2.0 monolayer 
(ML) (1.9 ML for the sample grown at 510°C) was then 
deposited with a calibrated rate of 0.1 ML/s at a con- 
trolled substrate temperature. A gradually changed InAs 
amount was achieved by stopping the substrate rotation. 
This method was widely used in studying the QD 
growth dynamics and to fabricate QD samples with low 
areal density [8,10,13,17]. The effective indium flux and 
real deposition amount could be calibrated based on the 
cosine law for certain configurations of the MBE source 
beam [18]. Growth interruption (GI) of 10 s was intro- 
duced after InAs layer deposition. A 100-nm GaAs cap- 
ping layer was then grown at 600°C. Details of the 
sample growth processes can be found in another study 
[16]. For further spectroscopy measurements, the sam- 
ples were cut into 16 pieces along the direction corre- 
sponding to which the InAs amount increased gradually. 



To evaluate the WL information, the relative reflectance 
difference in the sample surface plane, i.e., r/r = 
2(r [n0 ] -r [n _o])/(r [110 ] +r [n _ 0] ) , was measured with 
the RDS technique in ambient conditions. The setup of 
our RDS was reported elsewhere [19]. 

Results and discussion 

Figure la shows the intensity map of the second- 
derivative spectra obtained from the samples grown at 
530°C, in which the distinctive features of the GaAs 
band edge, LH and HH, of the WLs can be distin- 
guished. For the horizontal axis, the InAs deposition rate 
is calibrated with the cosine law mentioned above. One 
can see the LH and HH transition energies redshift al- 
most linearly up to an InAs deposition rate of 0.085 
ML/s (corresponding to an InAs amount of 1.7 ML), 
which is commonly known as the CT of InAs QD for- 
mation, then gradually saturate for further deposition. In 
order to have an intuitive understanding on the evolu- 
tion, we calibrate the WL thicknesses based on the tran- 
sition energies obtained from the spectra [20]. Figure lb 
gives the WL thickness evolution of the six samples. The 
WL thicknesses are signed with open symbols for 2D 
growth stage and solid symbols for 3D growth stage 
based on our previous results [16]. In general, two dis- 
tinct evolution processes can be discerned. An initially 
increased and then saturated evolution mode was 
observed in the samples grown at lower temperatures, 
while distinct concave-up features of the WL evolution 
for samples with relatively higher growth temperatures 
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Figure 1 Second-derivative spectra and WL thickness evolution, (a) The second-derivative spectra {d 2 p I d\ 2 ) of the series of samples grown 
at 530°C indicated by color contrast. The wavelengths of GaAs band edge, LH- and HH-related transitions in the WL are indicated by squares, 
circles, and triangles, respectively, (b) WL thickness evolution with InAs deposition amount for samples grown at different temperatures. Note that 
WL thicknesses are marked with open symbols at 2D growth stages, and solid symbols at 3D growth stage. The 2D to 3D transition points are 
determined in our previous work [16]. 
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appeared right after QD formation. They are unexpected 
based on thermodynamic understanding of the SK 
growth, for which a stable WL thickness is expected dur- 
ing and after QD formation. There is another decrease 
of the WL thickness for the last three pieces of each 
sample. The generation of large QDs, dislocation during 
the QD ripening process, and the enhanced indium ab- 
sorption abilities by them can be accounted for the de- 
crease [14,21]. 

In describing the WL growth dynamics, we consider 
three main InAs distribution processes. The newly 
deposited InAs can be incorporated into the WL and 
QDs or re-evaporated through the indium desorption 
process. Other processes, such as the formation of 
quasi-QDs or indium droplets, are neglected for their 
relatively lower rates [22,23]. In our previous works, the 
CTs of those samples were successfully determined, 
which enabled us to consider the 2D and 3D growth 
stages separately in this rate equation model [16]. For 
2D growth stage, the deposited material contributes to a 
pseudomorphic growth of the WL and the formation of 
QDs is neglected. Based on material balance, the rate 
equation can be written as follows: 



In previous works, for a WL thickness of 6, the instabil- 
ity of WL is commonly considered as (6 - 6 C ). The driv- 
ing force of QD growth, which is known as superstress', 
is defined as £ = (6 - 6 C ) I 0 C [27]. However, it is not 
suitable in describing our experimental results. The 
concave-up style of evolution shown in the upper panel 
of Figure lb means a non-zero QD growth rate when 
the WL thickness is slightly below the CT, or else, the 
WL thickness would not reduce below the CT in the 
presence of sufficient InAs supply. A non-stopping QD 
formation when the WL thickness is smaller than the 
CT is also documented in previous experiments [28-30]. 
Hence, the unstable part of WL is written as (6 - a6 c ), 
(0 < a < 1); correspondingly, the superstress is written as 
£ = (6 - a6 c ) I aO c (0 < a < 1). The QD formation and 
growth rate y is considered to be exponentially 
dependent on the superstress or y = b exp where b 
and were constant parameters [31,32]. Consequently, 
F QD can be written as follows: 



*2k T 

F QD = b(0 - a6 c ) — j— exp(-E di f/k B T)exp(fiZ). 



(4) 



dt 



,(t<t c ), 



(1) 



where 6 is the WL thickness, t c is the time of the 2D 
growth stage, and G is the InAs deposition rate. The in- 
dium desorption rate is presumed to be proportional to 
the InAs amount in WL, and r des represents the desorp- 
tion time constant. The indium desorption process is 
generally considered as thermally activated [23,24]; r des 
can be written as T des = ^ exp (j^j , where v 0 is a pre- 
exponential factor and E des is the activation energy of 
the indium desorption process. By solving Equation 1, 
the WL thickness versus growth time can be written as 
follows: 



6 = Gr des 1 - exp 



T-des 



.(t < t c ) 



(2) 



Above the CT, a large amount of QDs appears, and 
the newly deposited InAs are mainly consumed by them. 
The rate equation can be written as follows: 



dt 



F QD .(t>t c ) 



(3) 



Here F QD is used to represent the InAs consumption 
rate by the QD formation and growth [25]. F QD is deter- 
mined by the instability of WL and the material diffu- 
sion from WL to QDs [26]. The diffusion rate can be 
written as D In = (2k B T I h) exp {-E di f I k B T), [26] where 
k B is Boltzmann's constant, h is Planck's constant, T is 
the substrate temperature and E di f is the energy barrier. 



From Equation 3, it is clear that an equilibrium WL 
thickness is reached when the deposition rate equals to 
the WL consumption rate by QD formation and indium 
desorption. The last two processes show strong 
temperature dependence. Thus, in principle, WL growth 
has deposition rate and temperature dependences. WL 
growth dynamics at different conditions can be obtained 
from Equation 2 and by solving Equation 3 numerically. 
The calculation results for two different temperatures, 
and varied deposition rates are shown in Figure 2. The 
2D to 3D transition (where the WL thickness exceeds 
the CT) will not necessarily appear during the growth, 
depending on the deposition rate and time. Generally 
speaking, firstly, the WL thickness shows a nearly linear 
increase and then is saturated after QD formation. After 
that, the newly deposited InAs is mainly consumed by 
the formation of QDs. WL growth tends to reach equi- 
librium. The larger the InAs deposition rate is, the 
thicker the steady-state WL. Concerning the influence of 
growth temperature, distinct differences can be found by 
comparing the 3D growth stage in Figure 2a,b. For the 
low-temperature case, the WL thickness increases to 
equilibrium values which are always above the CT. On 
the other hand, for high- temperature grown samples, 
they show WL erosion at the beginning of the 3D 
growth stage, and WL thicknesses stable at values smal- 
ler than the CT. WL erosion disappears by increasing 
the deposition rates. According to Equation 3, the ap- 
pearance of WL erosion (d6 I dt < 0) resulted from the 
temperature dependence of QD formation rate, F QD . WL 
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Figure 2 Calculation results of WL growth dynamics for different InAs deposition rates at (a) 490°C and (b) 520°C. The horizontal and 
vertical dotted lines correspond to the critical thickness and t = 20 s, respectively. The insets of (a) and (b) zoom in the 3D growth stages. Here, 
we use b = 83 x 10" 9 , a = 0.85, /3 = 10, and E dif = 1.04 eV to solve the rate equations numerically. 



thickness would suffer a decrease if the deposition rate 
is not large enough to sustain the QD growth at the be- 
ginning of the 3D evolution stage. F QD then drops cor- 
respondingly according to its dependence on the 
superstress. It takes some time for the WL thickness to 
be stable to a certain value, for which the material de- 
position and the QD growth reach a balance. The bigger 
the deposition rate is, the thicker the stabilized WL. 
However, if the QD formation rate for the critical WL 
thickness is lower than the corresponding deposition 
rate, e.g., a lower growth temperature, WL thickness 
would keep on increasing after QD formation, which is 
the case in Figure 2a. 

To understand the experimental results, we have to 
apply the above-mentioned model to the non-rotating 
samples. The WL thicknesses after deposition can be 
known by calculating the resulting WL thicknesses with 
gradually changed InAs deposition rates and a given 



deposition time. Considering the WL evolution during 
GI, one can simply set G = 0 in Equations 2 and 3 and 
calculate with the WL thicknesses after deposition as 
initial values. The simulation results of WL evolution at 
different temperatures are shown in Figure 3a. One 
could see that the main features are well reproduced 
compared with the experimental results. The WL thick- 
ness shows a monotonic increase if the temperature is 
set at 490°C but shows a concave-up evolution for tem- 
peratures of 520°C and 540°C. According to the discus- 
sion above, we know that the slowed down increase 
observed at lower temperature is because of the depos- 
ition rate dependence of the equilibrium WL thickness. 
For higher growth temperatures, the elevated QD forma- 
tion rate at the beginning of the 3D growth stage led to 
WL erosion, which corresponds to the decrease of WL 
thickness on those non-rotating samples. WL thickness 
increases again with the deposition rate when the growth 
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Figure 3 Simulation of the WL evolution and fitting result of temperature dependence, (a) Simulation of WL evolution of the non-rotating 
samples at three different temperatures. The open and closed symbols stand for 2D and 3D growth stages, respectively. Here, we use 
b = 4 x 10~ 9 , a = 0.8, /3 = 10, and E dif = 1.04 eV to solve the rate equations numerically. The inset of (a) is the dependence of 3D growth time on 
InAs deposition rate for a given deposition time of 20 s. (b) The fitting result of the temperature dependence of WL thicknesses for the sample 
with a nominal InAs deposition amount of 1.45 ML. The inset shows the temperature dependence of desorption time constants for samples with 
different InAs deposition rates. 
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reaches equilibrium. It should also be noted that the 
simple equations do not reproduce the experimental 
results quantitatively because of their semi-empirical na- 
ture and the use of some adjustable parameters. 

We would like to comment on a special feature of 
those non-rotating samples. The material deposition rate 
changes gradually at different positions of a sample, 
which leads to the same behavior of deposition amount 
for a given growth time. Considering a weak dependence 
of the CT on deposition rate, one would expect that it 
takes different times at those positions of the sample to 
enter the 3D growth stage. The 2D growth time t 2 n can be 
calculated respectively from Equation 2 by taking 0 = 6 C . 
One then obtains the 3D growth time t 3D = t InAs - t 2 D- 
The inset of Figure 3a shows the 3D growth time with 
deposition rates. It should be noticed that at some posi- 
tions, they have very small values. Apparently, a near- 
zero 3D growth time cannot ensure an equilibrium 
quantum dot growth nor provide a steady-state WL 
thickness. It leads to stronger kinetic control characters 
on those samples. 

We come back to the 2D growth stage and to study 
another kinetic problem during growth - the indium de- 
sorption. In Figure lb, for samples with the same depos- 
ition amount but different growth temperatures, WL 
thicknesses varied a lot. From Equation 1, we know that 
it is because of the temperature dependence of indium 
desorption rate. From Equation 2 and from further con- 
sidering the WL evolution during GI, the resulting WL 
thickness in the 2D growth stage can be written as fol- 
lows: 

fl =Gr fc (l-exp(-^))exp(-£), (5) 

where t InAs is the InAs deposition time and t GI is the GI 
time. The kinetic parameter of indium desorption, E desf 
and v 0 , can be extracted from Equation 5 and Figure lb. 
We adopt the WL thicknesses of the first four pieces of 
each sample with effective InAs deposition amounts of 
1.14, 1.24, 1.34, and 1.45 ML to fit E des and v 0 , respect- 
ively. The obtained E des = 3.68 eV and v 0 are around 
5.5 x 10 22 . The activation energy is close to previously 
reported InAs decomposition energy and indium de- 
sorption activation energy from InGaAs [23,33]. We no- 
tice that the fitting v 0 is such a big number; v 0 stands for 
the attempt frequency of desorption, which is commonly 
known with the order of 10 12 to 10 14 s" 1 for desorption 
from metal and semiconductor surfaces. Such a big tran- 
sition frequency obtained here is also reported by other 
groups in investigating the InAs/GaAs QD desorption 
[33] or As desorption from GaAs surface [34]. It is con- 
sidered as physically achievable and could explain sev- 
eral characteristic features in InAs MBE growth [34]. 



The inset of Figure 3b shows the temperature depend- 
ence of the desorption life time {r des ) for samples with 
different InAs deposition amounts based on the fitting 
results. The time constants show a weak dependence on 
the indium flux but strongly decrease with increasing 
temperature; r des decreases from 1,063 s at 490°C to 35 s 
at 540°C for samples with a deposition amount of 
1.45 ML. The same strong dependence is also mentioned 
elsewhere [35]. Those time constants could be used to 
estimate the degree of desorption during the growth of 
InAs/GaAs(001) QDs at a certain temperature. 

Conclusions 

In conclusion, two kinds of WL evolution process of 
InAs/GaAs(001) QD system have been discussed based on 
RDS measurements and a rate equation model. They were 
well understood in the regime of material balance of WL 
growth/consumption and temperature dependence of QD 
formation. The concave-up style of evolution is also an 
evidence of a non-zero QD growth rate when the WL 
thickness was slightly lower than the critical value. We 
also gave a brief discussion on the indium desorption 
process during growth. Those results helped us in under- 
standing the kinetically controlled QD growth process. 
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